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Abstract 

A  model  equation  is  developed  to  predict  the  cell  voltage  versus  current  density  response  of  a  liquid  feed  direct  methanol  fuel  cell 
(DMFC).  The  equation  is  based  on  a  semi-empirical  approach  in  which  methanol  oxidation  and  oxygen  reduction  kinetics  are  combined 
with  effective  mass  transport  coefficients  for  the  fuel  cell  electrodes.  The  model  equation  is  validated  against  experimental  data  for  a 
small-scale  fuel  cell  and  is  applicable  over  a  wide  range  of  methanol  concentration  and  temperatures. 
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1.  Introduction 

Characterisation  of  fuel  cells  frequently  uses  large  and 
complex  computer  models,  based  on  minute  details  of  cell 
component  design  (physical  dimensions,  materials,  etc.) 
along  with  chemical  and  physical  considerations  (transport 
phenomena,  electrochemical  kinetics,  etc.).  The  codes,  of¬ 
ten  proprietary,  needed  in  the  design  and  development  of 
fuel  cells,  are  cumbersome  and  time-consuming  for  use  in 
system  analysis  models.  Simpler  approaches  are  normally 
used  for  system  studies.  Another  approach,  which  is  not 
time  and  cost  efficient,  would  be  to  conduct  appropriate  tests 
at  every  condition  expected  to  be  analysed  in  the  system. 
Alternatively,  it  is  prudent  to  develop  correlations,  based  on 
thermodynamic  modelling,  that  describe  cell  performance 
as  operating  conditions,  such  as  temperature  and  pressure 
are  changed  [1].  Thermodynamic  modelling  is  used  to  de¬ 
pict  the  equations  so  that  only  a  limited  number  of  tests  are 
needed  to  define  design  constants  within  the  equation. 

In  the  development  of  model  equations  to  predict  per¬ 
formance  of  polymer  electrolyte  fuel  cells  a  number  of  ap¬ 
proaches  using  empirical  models  have  been  attempted  [2-7] . 
In  many  cases  an  excellent  agreement  between  the  model 
and  experimental  data  is  achieved  by  adjusting  appropriate 
coefficients/parameters  in  the  model  equations.  However,  a 
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major  limitation  is  that  the  coefficients  do  not  follow  any 
specific  trends  with  fuel  cell  operating  variables  and  do  not 
allow  any  physically  real  interpretation  of  the  model.  Such 
an  approach  has  also  been  applied  to  the  case  of  the  direct 
methanol  fuel  cell  (DMFC)  [8]. 

Direct  methanol  fuel  cells  have  reached  a  significant 
stage  of  development  with  high  power  density  performance 
now  achieved.  Probably  without  exception  the  research 
and  development  of  the  DMFC  to  reach  high  performance 
targets  has  proceeded  with  anodes  of  the  binary  alloy  of 
Pt  and  Ru.  In  general,  the  performance  of  the  DMFC  de¬ 
pends  on  the  MEA  construction  method,  the  materials  and 
the  operating  conditions  in  the  cell.  It  is  well  documented 
that  increasing  cell  temperature  above  approximately  60  °C 
causes  a  significant  increase  in  power  performance  and  that 
at  90  °C  and  beyond  power  densities  above  200mWcm~2, 
are  achievable.  This  level  of  power  performance  is  typi¬ 
cally  achieved  using  a  pressurised  cathode  oxygen  or  air 
(2-5  bar).  The  use  of  oxygen  gives  the  best  performance  but 
most  terrestrial  applications  require  the  use  of  air.  In  addi¬ 
tion,  many  potential  applications  of  the  DMFC  will  require 
air  under  almost  ambient  conditions  of  temperature  and 
pressure. 

In  this  paper,  we  develop  a  semi-empirical  model  for  the 
cell  voltage,  current  density  response  of  the  DMFC  which  is 
based  on  Tafel  type  kinetics  for  methanol  oxidation  and  oxy¬ 
gen  reduction  and  on  measured  mass  transport  coefficients. 
The  model  equation  is  validated  against  experimental  data 
obtained  for  a  small-scale  fuel  cell. 
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Fig.  1.  DMFC  small-scale  cell  flow  bed  design  and  manifold  arrangement. 


2.  Experimental 

Tests  on  the  DMFC  were  performed  with  a  cell  with 
a  cross-sectional  area  of  9  cm2.  The  cell  consisted  of  two 
non-porous  graphite  blocks  with  a  series  of  parallel  chan¬ 
nels  machined  into  the  block  for  the  flow  of  methanol  and 
oxygen/air  (Fig.  1). 

The  cell  was  fitted  with  a  membrane  electrode  assem¬ 
bly  (ME A)  sandwiched  between  the  two  graphite  blocks. 
The  cell  was  held  together  between  two  plastic  insulation 
sheets  and  two  stainless  steel  backing  plates  using  a  set 
of  retaining  bolts  positioned  around  the  periphery  of  the 
cell.  Electrical  heaters,  supplied  by  Watson  Marlow,  were 
placed  behind  each  of  the  graphite  blocks  to  heat  the  cell 
to  the  desired  operating  temperature.  The  graphite  blocks 
were  also  provided  with  electrical  contacts  and  small  holes 
to  accommodate  thermocouples.  The  fuel  cells  were  used 
in  a  simple  flow  rig,  which  consisted  of  a  Watson  Marlow 
peristaltic  pump  to  supply  aqueous  methanol  solution,  from 
a  reservoir,  to  a  Eurotherm  temperature  controller  to  main¬ 
tain  the  cell  at  a  constant  temperature.  Air  was  supplied 
from  cylinders,  at  ambient  temperature,  and  the  pressure 
regulated  by  needle  valves. 

MEAs  studied  in  this  work  were  made  in  the  following 
manner:  the  anode  consisted  of  a  Teflonised  (20%)  car¬ 
bon  cloth  support  (E-Tek,  type  ‘A’),  of  0.3  mm  thickness, 
upon  which  was  spread  a  thin  (diffusion  layer)  layer  of 
uncatalysed  (ketjenblack  600)  10wt.%  teflonised  carbon. 
The  catalysed  layer,  unless  otherwise  specified  consisted 
of  35wt.%  Pt-15wt.%  Ru  (2  mg  cm-2  metal  loading)  dis¬ 


persed  on  carbon  (ketjen)  and  bound  with  bound  with 
10wt.%  Nation®  from  a  solution  of  5wt.%  Nation®  dis¬ 
solved  in  a  mixture  of  water  and  lower  aliphatic  alcohol’s 
(Aldrich),  was  spread  on  this  diffusion  backing  layer.  A 
thin  layer  of  Nation®  solution  was  spread  onto  the  sur¬ 
face  of  each  electrode.  Electrode  preparation  are  described 
in  detail  elsewhere  [10].  The  cathode  was  constructed 
using  a  similar  method  as  for  the  anode,  using  a  thin  dif¬ 
fusion  layer  bound  with  10wt.%  PTFE,  and  lmgcm-2 
Pt  black  with  10wt.%  Nation®  in  the  catalyst  layer.  The 
electrodes  were  placed  on  opposite  sides  of  a  pre-treated 
Nation®  117  membrane.  This  pre-treatment  involved  boil¬ 
ing  the  membrane  for  1  h  in  5  vol.%  H2O2  and  1  h  in 
1  mol  dm-3  H2SO4  before  washing  in  boiling  Millipore 
water  (>18m£2)  for  2h  with  regular  changes  of  water. 
The  assembly  was  hot-pressed  at  100  kg  cm-2  for  3  min 
at  135  °C. 

Cell  voltage  versus  current  density  response  was  mea¬ 
sured  galvanostatically,  by  incrementally  increasing  the 
current  from  open  circuit  and  measuring  the  cell  potential 
and  then  reducing  the  current  incrementally  again  measur¬ 
ing  the  cell  voltage.  The  data  reported  here  were  obtained 
with  one  MEA,  unless  otherwise  stated.  The  MEA  was  con¬ 
ditioned  before  use  in  two  stages:  for  48  h  in  the  test  cell, 
in  2.0  mol  dm-3  methanol  solution  at  75  °C,  and  then  by 
maintaining  the  cell  with  an  applied  load  of  100  mA  cm-2 
for  several  hours.  This  pre- treatment  resulted  in  stable  per¬ 
formance  under  continuous  operation.  Data  used  in  this 
empirical  modelling  has  been  presented  earlier  in  a  previous 
publication  [11]. 
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3.  Overall  cell  voltage 

It  is  well  known  that  useful  work  (electrical  energy)  is 
obtained  from  a  fuel  cell  only  when  a  reasonable  current 
is  drawn,  and  that  the  actual  cell  potential  falls  below  its 
equilibrium  potential  because  of  several  irreversible  losses: 
catalyst  polarisation,  internal  resistance  and  mass  transport. 

In  general,  the  overall  cell  voltage  ECQ\\  has  several  com¬ 
ponents: 

(i)  Equilibrium  potentials,  determined  from  the  Nernst 
equation,  of  cathode  and  anode  reactions  (E^  and  E ^), 
which  make  up  the  reversible  cell  voltage,  EQ. 

(ii)  Overpotentials  at  the  cathode  and  the  anode,  r /a  and  ijc- 

(iii)  Ohmic  losses  in  the  electrolytes  (/^electrolyte),  the  cell 
separator  (/7?sep),  the  electrodes  and  in  the  connections 
from  the  power  supply  to  the  electrodes  (/^rr/^rJ. 

From  the  above,  the  cell  voltage  can  be  presented  in  the 
form: 


^cell  =  -  Et  -  l^cl  -  M  -  IRc irc  “  -^dre 

/^electrolyte  /^separator  (1) 

The  overpotentials  at  electrodes  arise  from  various  polar¬ 
isation  phenomena  and  increase  in  magnitude  as  the  rate  of 
reaction,  or  the  current  density,  increases.  The  electrode  de¬ 
sign  and  material  should  also  facilitate  efficient  gas  release 
from  its  surface  to  ensure  that  the  IR  drop  in  the  electrolyte 
and  the  “bubble  polarisation”  is  low.  The  material  of  the  cur¬ 
rent  connectors  should  be  highly  conducting  and  designed 
to  the  thickness  required  to  carry  anticipated  maximum  cell 
currents.  The  electrode  material  should  also,  if  possible  be 
highly  conducting  and  be  of  the  required  thickness.  A  low 
electrical  resistance  of  the  electrodes  is  also  desirable  to 
achieve  a  uniform  potential  distribution  over  the  face  of  the 
electrode.  The  practical  size  of  the  electrodes  is  often  lim¬ 
ited  by  these  two  considerations. 

3.1.  Electrode  polarisation 


The  objective  of  a  fuel  cell  is  to  convert  as  much  of  the 
fuel  as  possible  to  the  product,  the  maximum  amount  of 
the  work  which  may  be  derived  from  the  fuel  will  be  that 
corresponding  to  —  AG  after  conversion,  i.e.  at  the  product 
exit  of  the  fuel  cell.  There,  —  AG  will  be  less  negative  than 
the  value  at  reactant  inlet,  since  the  concentration  of  the 
reactant  will  have  fallen  and  correspondingly  that  of  the 
product  will  have  risen.  Thus,  the  maximum  voltage  Em ax 
that  a  fuel  cell  can  develop  under  reversible  conditions  at 
real  conversions  is  given  by 


11 F Em3X 


-AG  =  -AG°t  +/?rin 


where  n  is  the  number  of  electrons  exchanged  per  molecule 
of  reactant,  F  the  Faraday  constant,  R  the  gas  constant, 
—  A  GP  is  the  Gibbs  energy  of  the  reaction  under  standard 


concentration  (activity)  conditions  at  operating  temperature 
T  (K);  Cr  and  Cp  are  the  ratios  of  concentrations  (or  more 
exactly  activities)  of  reactant  and  product  leaving  the  fuel 
cell  to  those  in  the  standard  state,  and  m  is  the  reaction  order 
of  product  relative  to  that  of  the  reactant. 

The  polarisation  losses  are  determined  by  the  kinetics  of 
the  electrode  reactions,  by  the  physical  structure  (geometry) 
of  the  cell  and  by  the  type  of  electrode  used.  From  a  practi¬ 
cal  point  of  view,  the  “overvoltage”  is  the  voltage  difference, 
which  is  measured  between  the  “open  circuit  voltage”  and 
the  “terminal  voltage”  under  the  conditions  of  current  flow¬ 
ing  in  either  direction.  Tafel’s  equation  is  typically  used  to 
represent  overvoltage,  rj: 

r)  =  a  +  b  In  jc  =  a  +  b  In  (nFk)  (3) 


where  jc  is  the  cathodic  current  density  (A  cm-2),  k  the 
rate  expressed  (mol  cm-2  s-1),  n  the  number  of  electrons 
successively  transferred  per  reacting  molecule,  and  a  and  b 
are  constants  for  a  given  reaction  and  substrate,  at  a  given 
temperature.  The  Tafel  constant  b  is 


RT 

b  = - 

olF 


where  ac  is  the  cathodic  transfer  coefficient. 

The  negative  sign  indicates  that  the  cathodic  current  den¬ 
sity  increases  as  V  decreases.  In  the  anodic  direction  (e.g. 
hydrogen  oxidation),  the  corresponding  equation  has  con¬ 
stant  ot!  and  b\  where  b'  is  equal  to  RT/a'F.  For  cases  com- 


monly  encountered: 

0ic  =  n  /3(/trds) 

(5) 

aa  =  n"  +  (1  -  P)(nvds) 

(6) 

where  aa  is  the  corresponding  transfer  coefficient,  j3  is 
called  the  symmetry  factor,  n'  the  number  of  electrons 
which  are  transferred  in  the  reaction  sequence  before  each 
single  rate-determining  step  (rds)  in  the  cathodic  direction, 
in  which  (nrds)  electrons  are  transferred,  and  n"  is  the  num¬ 
ber  similarly  transferred  in  the  anodic  direction;  (nrds)  can 
only  be  1  (for  an  electrochemical  rds)  or  0  (for  a  chemical 
rds).  Thus 


ac+aa  =  n  (7) 

where  n  is  the  number  of  electrons  transferred  in  the  overall 
process  per  unit  rds. 

The  net  electrode  reaction  is  the  sum  of  the  anodic  (pos¬ 
itive  current)  and  cathodic  (negative  current)  processes  and 
is  given  by  the  Butler- Volmer  expression: 

j  =  ;0[e(a *Frl/RT)  -  Q-GcFri/RT)}  (g) 

3.2.  Empirical  equations  for  PEMFC 

Srinivasan  et  al.  [2]  showed  that  it  is  possible  to  use  a 
simple  equation  to  describe  the  cell  voltage  versus  current 
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density  behaviour  for  PEMFCs  in  the  activation  and  Ohmic 
controlled  current  density  region: 

E  =  E0-b\ogj-RJ  (9) 

E0  =  Ev  +  b  log  jo  (10) 

where  Ev  is  the  reversible  cell  potential,  b  is  the  Tafel  slope 
for  oxygen  reduction  and  R  is  the  Ohmic  resistance  of  the 
solid  polymer  electrolyte. 

Using  Eq.  (1)  it  was  shown  that  as  current  density  is 
increased  the  predicted  cell  potential,  with  the  appropriate 
coefficients,  decreases  much  less  rapidly  than  observed  [3]. 
To  increase  the  reliability  of  the  aforementioned  equation 
Kim  et  al.  introduced  a  more  complex  equation: 

E  =  E0  -  b  log  i  -  Rti  -  mtni  (1 1) 

where  m  and  n  are  parameters  that  account  for  the  “mass 
transport  overpotential”  as  a  function  of  current  density. 

Parameter  m  affects  both  the  slope  of  the  linear  region  of 
the  E  versus  j  plot  and  the  current  density  at  which  there  is 
a  departure  from  the  linearity.  The  value  of  n  has  a  predom¬ 
inant  effect  in  the  mass  transport  limitation  region.  Kim  re¬ 
ported  that  despite  the  excellent  fit  when  applying  Eq.  (11) 
to  experimental  data  they  were  not  able  to  find  a  theoreti¬ 
cal  explanation  for  the  last  term  of  their  equation  [3].  Lee 
and  Lalk  proposed  [4-6]  a  slightly  modified  form  of  Kim’s 
equation  by  adding  the  ratio  of  cell  pressure  to  oxygen  par¬ 
tial  pressure 

V  =  E0  —  b  log  I  —  RqI  —  m  exp,?/  —  b  log  (12) 

More  recently  Squadrito  et  al.  [7]  used  Eq.  (1 1)  as  a  start¬ 
ing  point  and  he  tried  to  analyse  the  different  contributions 
to  the  mass  transport  limitation.  He  considered  two  differ¬ 
ent  contributions;  an  Ohmic  one  which  is  independent  of 
the  current  density  and  a  non- Ohmic  one  which  varies  with 
current  density  and  produced  an  equation  in  the  form: 

E  =  E0  -  b  log  i  -  Rei  +  aik  ln(l  -  (ii)  (13) 

where  a ,  k  and  /3  are  coefficients. 

The  logarithmic  term  ln(l  —  /3i)  introduces  a  limit  to  the 
available  current  density.  For  k  =  1,  a  has  the  same  dimen¬ 
sion  of  Rq  and  can  be  interpreted  as  an  additional  resistance 
term  due  to  the  overall  mass  transport  limitation.  The  pa¬ 
rameter  k  mainly  influences  the  point  at  which  there  is  a 
departure  from  linearity  and  a  determines  the  shape  of  the 
curve  at  high  current  densities. 

Both  empirical  Eqs.  (12)  and  (13)  include  several  coef¬ 
ficients,  which  we  have  used  to  give  good  agreement  with 
DMFC  fuel  cell  polarisation  data  [8].  However,  as  with  the 
equivalent  PEM  cell  model,  changes  in  the  operating  con¬ 
ditions,  e.g.  methanol  concentration,  produced  different  val¬ 
ues  of  the  empirical  coefficients  and  the  predictability  of  the 
model,  outside  a  narrow  parameter  domain,  was  lost.  Thus,  it 
is  important  to  establish  models  that  incorporate  the  physics 
and  chemistry  of  the  system  to  enable  good  prediction  of 
behaviour  over  a  broad  range  of  parameters  and  variables. 


4.  An  empirical  equation  for  modelling  liquid  DMFC 
behaviour 


The  DMFC  is  a  complex  porous  electrode  system  due 
to  the  thin  three-dimensional  structure  of  the  electrocatalyst 
regions,  the  porous  diffusion  layers  and  the  generation  of 
carbon  dioxide  gas.  The  DMFC  exhibits  limiting  current  be¬ 
haviour  due  to  mass  transport  limitations  of  methanol  supply 
to  the  anode  catalyst  caused  by  the  following. 

(i)  Hydrodynamic  mass  transfer  from  the  feed  flow  to  the 
surface  of  the  surface  of  the  carbon  cloth  backing  layer, 
where  the  gas  bubbles  are  released  into  the  flowing 
methanol  solution. 

(ii)  Diffusion  mass  transfer  in  the  carbon  cloth  and  anode 
carbon  diffusion  layer. 

(iii)  Diffusion  mass  transfer  in  the  catalyst  layer. 


Mass  transport  limiting  currents,  place  a  potentially  large 
restriction  on  the  performance  of  the  DMFC  and  indicates 
the  need  for  effective  electrode  design  and  suitable  mass 
transport  analysis.  The  combined  effect  of  (i)  and  (ii)  is 
simply  expressed  as 


1  _  i  i 

kQ  ff  ki  kc  i 


(14) 


where  the  effective  mass  transfer  coefficients,  keff,  can  be 
determined  at  the  limiting  current  from 


j 

uFCme 


(15) 


where  Cme  is  the  methanol  concentration. 

ki  is  the  hydrodynamic  mass  transport  coefficient,  given 
by 

h  =  r(\j)P  (16) 


where  p  =  0.32  and  r  =  1.87  x  10-6  for  the  fuel  cell 
considered  in  this  paper  [11]. 

The  value  of  ki  depends  on  the  hydrodynamics  in  the 
channel  as  well  as  on  gas  bubble  liberation  at  the  surface. 
kc i  is  the  carbon  cloth  mass  transport  coefficient,  given  by 

kci  =  k°cle™  (17) 


where  ec  is  the  liquid  voltage,  m  =  1.5,  and  k^  =  (Dmeo h/ 
lc i),  where  lc\  is  the  effective  backing  layer  thickness. 

In  practice,  mass  transport  in  the  backing  layers  is  unlikely 
to  be  governed  solely  by  diffusion  as  liquid  and  gas  flow 
in  the  channels  and  gas  evolution  will  cause  “mixing”  in 
regions  facing  the  flow  channel  and  thus  reduce  the  effective 
thickness  of  the  carbon  backing. 

The  porous  electrocatalyst  layer  is  initially  considered  as 
a  thin  pseudo  two-dimensional  structure  where,  due  to  the 
potential  distribution,  highest  activity  is  at  a  region  next  to 
the  membrane  surface. 

The  above  model  Eq.  (15)  predicts  that  mass  transfer  co¬ 
efficients  for  the  DMFC  will  exhibit  a  broad  maximum  in 
value  with  increasing  methanol  concentration.  The  model 
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Fig.  2.  Effective  mass  transfer  coefficients  for  the  direct  methanol  fuel  cell  determined  at  limiting  current  conditions.  Points  are  experimental  data. 


gives  reasonable  agreement  with  experimentally  measured 
mass  transport  coefficients  for  the  DMFC  (Fig.  2).  The 
mass  transport  coefficients  are  in  the  range  of  approximately 
(2. 5-6.0)  x  10  6 ms  1  for  a  MEA  under  consideration  in 
this  paper. 

As  a  model  for  methanol  oxidation,  Tafel  type  kinetics  is 
chosen: 


(18) 


where  j0  is  the  exchange  current  at  the  reference  concentra¬ 
tion,  and  a  the  transfer  coefficient,  a  refers  to  the  condition 
at  the  anode  catalyst  surface  and  N  is  an  order  of  reaction. 

Recent  data  obtained  in  this  laboratory  has  confirmed  the 
applicability  of  the  Tafel  approximation  for  methanol  oxi¬ 
dation,  for  a  specific  electrode  assemble  considered  in  this 
paper  [12].  Although  we  use  the  Tafel  approximation  for 
electrode  kinetics  we  develop  an  empirical  equation  for  the 
full  range  of  current  densities;  from  zero  to  high  values.  This 
clearly  is  an  approximation  at  low  current  densities  (over¬ 
potential)  but  this  region  of  operation  is  not  of  real  practical 
interest  in  fuel  cell  operation. 


We  represent  mass  transport  in  the  anode  side  of  the  cell 
using  the  effective  mass  transport  coefficient 

j  =  keffnF(Cu E  -  CJra)  (19) 


By  re-arranging  the  previous  equation  the  overpotential 
becomes 
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By  combining  Eqs.  (19)  and  (20)  we  obtain 
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(21) 


Using  Eq.  (21)  for  the  anode  we  can  obtain  an  expression 
for  the  cell  voltage  by  introducing  the  kinetic  expression  for 
oxygen  reduction: 

J  —  Joe 

where  po  is  the  partial  pressure  for  oxygen. 
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Fig.  3.  Comparison  between  experimental  data  [11]  and  empirical  equation  based  prediction  for  a  cell  operated  with  0.125  M  methanol  solution  supplied 
at  a  rate  of  1.12cm3  min-1  with  air  fed  cathodes  pressurised  at  2bar  (cell  temperatures:  (O)  343.15;  (♦)  348.15  K;  (A)  353.15  K;  (■)  358.15  K;  (O) 
363.15  K). 

The  cathode  overpotential  is  then 


(23) 


0.8 


0.7 


0.6 


0.5 

> 

D) 

B  0.4 

o 

> 


0.2 


0.1 


0 

0  0.005  0.01  0.015  0.02  0.025  0.03  0.035  0.04  0.045 

Current  Density/  A  cm'2 

Fig.  4.  Comparison  between  experimental  data  [11]  and  empirical  equation  based  prediction  for  a  cell  operated  with  0.25  M  methanol  solution  supplied 
at  a  rate  of  1.12cm3  min-1  with  air  fed  cathodes  pressurised  at  2bar  (cell  temperatures:  (♦)  343.15;  (A)  348.15  K;  (•)  353.15  K;  (O)  358.15  K;  (O) 
363.15  K). 
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where  k\o  is  the  mass  transport  coefficient  for  the  cathode 
side  of  the  cell. 

Combining  Eqs.  (21)  and  (23),  we  obtain  an  expression 
for  the  cell  voltage,  /sceii,  as 
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The  above  equation  is  a  similar  form  to  the  Squadrito 
equation  with  k  =  0. 


5.  Model  validation 
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In  many  cases,  the  reduction  of  oxygen  will  not  proceed 
under  mass  transport  limitations  and  we  can  write  Eq.  (24) 
in  a  simplified  form  as 

£ceii  =  Eo~  ^ceii  log  j  -  RQj  +  C\  ln(l  -  C2j ) 
where 
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Eq.  (25)  above  was  used  to  model  the  cell  voltage  response 
of  a  small-scale  single  DMFC  cell  [11].  Model  predictions 
are  presented  for  four  different  aqueous  methanol  solution 
concentrations  namely  0.125,  0.25,  0.5  and  0.75  M  and  for  a 
range  of  cell  operating  temperatures.  As  it  can  be  seen  from 
Figs.  3-6  that  the  predictions  are  good.  The  advantage  of 
the  model  equation  lies  in  the  ability  to  follow  exactly  the 
voltage  profile  in  the  limiting  current  region. 

A  few  comments  are  necessary  regarding  the  limiting 
current  operating  region  of  the  DMFC.  As  can  be  seen 
from  the  figures  there  are  two  types  of  variation  in  voltage 
with  current  density,  depending  on  the  cell  operating  con¬ 
ditions.  One  is  a  steep  fall  in  voltage  with  current  density 
and  is  typical  of  cells  operated  under  severe  mass  trans¬ 
port  limitations  or  under  poor  reaction  kinetics  (low  tem¬ 
peratures,  low  methanol  concentrations,  or  combination  of 
both).  The  other  pattern  is  a  more  gradual  fall  in  voltage 
and  is  similar  to  that  commonly  seen  for  hydrogen  fuel 
cells. 


0  0.01  0.02  0.03  0.04  0.05  0.06  0.07  0.08  0.09  0.1 


Current  Density/  A  cm'2 


Fig.  5.  Comparison  between  experimental  data  [11]  and  empirical  equation  based  prediction  for  a  cell  operated  with  0.5  M  methanol  solution  supplied 
at  a  rate  of  1.12cm3  min-1  with  air  fed  cathodes  pressurised  at  2bar  (cell  temperatures:  (♦)  343.15;  (•)  348.15  K;  (A)  353.15 K;  (O)  358.15  K;  (O) 
363.15  K). 
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Fig.  6.  Comparison  between  experimental  data  [11]  and  empirical  equation  based  prediction  for  a  cell  operated  with  0.75  M  methanol  solution  supplied 
at  a  rate  of  1.12cm3  min-1  with  air  fed  cathodes  pressurised  at  2bar  (cell  temperatures:  (♦)  303.15;  (A)  313.151  K;  (O)  323.15  K;  (A)  333.15  K;  (O) 
343.15  K). 
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Fig.  7.  The  effect  of  varying  parameter  C2  of  Eq.  (21)  in  the  predicted  cell  response  profiles.  C2  varies  between  4.8  and  10.3  cm2  A  1  ( E0  =  0.3  V, 
b  =  0.12 Vdec-1,  R  =  0.0001  f2cm2,  Cx  =  0.059V). 
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A  common  weakness  of  many  models  for  the  DMFC  is 
that  they  are  only  valid  for  narrow  operating  conditions.  The 
current  equation  is  applicable  over  a  wide  range  of  operat¬ 
ing  conditions.  This  is  illustrated  in  Fig.  7  where  we  have 
used  a  set  of  coefficients  representing  one  of  our  data  sets 
and  have  varied  parameter  C2  of  Eq.  (21)  with  all  other 
parameters  remaining  constant.  This  data  shows  the  ability 
of  the  present  equation  to  predict  the  cell  potential,  current 
density  profiles.  In  principle,  this  means  that  it  can  be  used 


for  a  wide  range  of  membrane  electrode  assembly  materials 
and  cell  configurations  that  affect  mass  transfer;  for  exam¬ 
ple  using  different  backing  layers,  carbon  powders,  etc.  [8]. 
Alternatively,  this  data  of  Fig.  7  demonstrates  the  effect  of 
varying  the  methanol  concentration,  leaving  all  other  oper¬ 
ating  parameters  constant  [8,9]. 

Fig.  8  demonstrates  the  effect  of  decreasing  the  cell  re¬ 
sistance  in  Eq.  (21),  which  would  be  achieved  by  chang¬ 
ing  the  membrane  thickness  or  using  alternative  membrane 
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Fig.  8.  The  effect  of  decreasing  cell  resistance  parameter  Re  on  cell  voltage  R  varies  from  0.0001  to  2.2001  Q  cm2  (E0  =  0.3  V,  b  =  0.12  V  dec-1 
Ci  =  0.059  V,  C2  =  8  cm2  A-1). 
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materials  [10].  As  expected,  as  resistance  increases,  the  mass 
transfer  limitations  remain  essentially  unaffected,  but  the 
voltage  declines  more  rapidly. 

Fig.  9  demonstrates  the  effect  of  decreasing  parameter  Cy 
with  all  other  parameters  remaining  constant.  This  situation 
represents,  in  general,  limitations  related  to  reaction  kinetics 
which  are  attributable  to  a  variety  of  reasons,  e.g.  decrease 
in  cell  temperature  or  poorer  catalyst  performance. 

C2  varies  from  0.059  to  0.334  cm2  A-  (E0  =  0.3  V,  b  = 
0.12V dec-1,  R  =  0.0001  ^cm2  Cx  =  0.059V,  C2  = 
8  cm2  A”1). 

6.  Conclusions 

A  semi-empirical  equation,  based  on  Tafel  type  kinetics 
and  measured  mass  transport  coefficients,  has  been  devel¬ 
oped  to  predict  the  voltage  response  of  liquid  fed  direct 
methanol  fuel  cells  over  a  wide  range  of  operating  condi¬ 
tions.  The  proposed  equation  is  valid  even  in  the  case  of  very 
low  current  densities  caused  by,  for  example,  the  use  of  di¬ 
lute  methanol  solutions  or  low  cell  temperatures.  A  further 
paper  uses  the  model  to  determine  appropriate  kinetic  and 
resistance  parameters  for  the  DMFC  using  data  from  two 
ME  A  structures. 
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